Spin- Vacancy-Induced Long-Range Order in a New Haldane-Gap Antiferromagnet 
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Magnetic susceptibility, high-field magnetization and inelastic neutron scattering experiments 
are used to study the magnetic properties of a new S = 1 quasi- 1-dimensional antiferromagnet 
PbNi2V20g. Inter-chain interactions are shown to be almost, but not quite, strong enough to destroy 
the nonmagnetic singlet ground state and the energy gap in the magnetic excitation spectrum. 
Substituting nonmagnetic Mg 2+ (S = 0) ions for Ni 2+ (S = 1) induces a magnetically ordered state 
at low temperatures. To our knowledge this is the first observation of doping-induced long-range 
order in a Haldane-gap system. 

PACS numbers: 75.30.Kz, 75.30. Hx, 75.40. Cx, 75.40.Gb 



Unlike its classical counterpart, the quantum Heisen- 
berg antiferromagnet (HAF) has a nonmagnetic spin- 
singlet ground state in one-dimensional (ID) space. A 
pioneering work by Haldane |jj led to a breakthrough in 
our understanding of this remarkable phenomenon, and 
of low-dimensional magnetism in general. Haldane sug- 
gested that the actual ground states are totally differ- 
ent for integer- and half-integer-spin systems. The half- 
integer ID HAF shows quasi-long-range order (power-law 
spin correlations). True long-range order can be induced 
by the slightest perturbations. In contrast, for integer 
spins the ground state singlet is robust. Intrachain spin 
correlations decay exponentially (true short-range order) 
and the excitation spectrum has a nonzero energy gap. 
The system remains nonmagnetic even in the presence of 
weak magnetic anisotropy or residual 3-dimensional (3D) 
interactions. The spectacular properties of the integer- 
spin quantum HAF have been extensively studied theo- 
retically and experimentally, and are by now rather well 
understood. The current issue in quantum magnetism 
is how the singlet ground state is finally destroyed and 
long-range order (LRO) is formed when perturbations to 
the idealized Hamiltonian become sufficiently strong. 

Among the factors that favor an ordered ground state 
are single-ion anisotropy and 3D interchain coupling in 
quasi-lD HAF systems. Increasing either of these pa- 
rameters beyond some critical level leads to a quantum 
phase transition to a state with a Neel-like order param- 
eter. This effect has been extensively studied and the 
corresponding phase diagram is known from numerical 
simulations As anticipated by Shender and Kivel- 
son ||, another potential source of magnetic LRO is the 
introduction of spin- vacancies in the quantum chains. In- 



deed, replacing some of the S — 1 ions by S — sites 
severs the chains and liberates S = 1/2 end-chain de- 
grees of freedom Under certain circumstances these 
free spins can order in 3D space. To date however this 
phenomenon has not been observed in any material, and 
is poorly understood theoretically. In the present paper 
we report the first observation of spin-vacancy-induced 
ordering in a new Haldane-gap compound PbNi 2 V 2 8 , 
that sports a peculiar helical arrangement of magnetic 
ions and substantial (yet subcritical) interchain interac- 
tions. 

The exact crystal structure of PbNi 2 V 2 8 has not been 
reported yet [js]]. Preliminary X-ray diffraction measure- 
ments on powder samples prepared by conventional solid- 
state reaction suggest that this material is isomorphous 
with SrNi 2 V 2 8 || (Fig. 0, inset). SrNi 2 V 2 8 is tetrag- 
onal (space group 14 i cd) with lattice parameters a 
12. If) 1 7 A and c = 8.3247 A. We found that PbNi 2 V 2 8 
has the same crystal symmetry and very similar lattice 
parameters: a = 12.249(3) A and c = 8.354(2) A. In the 
two isostructural compounds slightly distorted NiOe oc- 
tahedra are edge-shared around the four-fold screw axis 
along the c direction. All the S = 1 Ni 2+ ions, as well as 
all the nearest neighbor (NN) Ni 2+ -Ni 2+ bonds are equiv- 
alent. The screw-chains are separated by VO4 tetrahedra 
and Sr 2+ or Pb 2+ ions. Intrachain nearest-neighbor Ni- 
Ni antiferromagnetic interactions are expected to domi- 
nate. The V 5+ sites are presumed nonmagnetic. 

To characterize the new Pb-based system, bulk mag- 
netic measurements were performed on highly aligned 
polycrystalline samples. These were prepared by apply- 
ing a magnetic field of 9 T at room temperature to a 
mixture of PbNi 2 V 2 8 powder and epoxy resin (Stycast 



*E-mail:uchinokura@ap.t.u-tokyo. ac.jp 



1 



1266). The Ni-chains tend to align themselves along the 
applied magnetic field. A c-axis alignment of better than 
3° (FWHM) was achieved, as determined from measuring 
the rocking curve of the 004 Bragg reflection. Magnetic 
susceptibility was measured with a commercial SQUID 
magnetometer. The experimental x(T) curves, shown in 
Fig. [j], have a broad maximum around 120 K for for both 
H _L c and H \\ c. Below this temperature, x(T) de- 
crease exponentially to almost zero. The residual low-T 
susceptibility can be attributed to the presence of para- 
magnetic impurities. A straightforward analysis of the 
data collected in the range < T < 40 K yields the 
thermal activation energies: A = 29.4 K = 2.53 meV for 
H || c, and A = 27.8 K = 2.39 meV for a transverse 
field, respectively. In addition to the susceptibility mea- 
surements, high-field magnetization data were taken on 
an aligned polycrystalline sample in magnetic fields of up 
to 40 T, using the pulsed magnetic facility at ISSP. Fig- 
ure^ shows magnetization curves measured at T = 4.2 K. 
At high fields an abrupt change of slope in M(H) is ob- 
served for both H \\ c and H _L c. The corresponding 
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FIG. 1. Temperature dependence of magnetic susceptibil- 
ity measured in PbNi2V20s. The solid line shows data for 
a nonaligned powder sample. The symbols are for highly 
aligned polycrystalline samples. The magnitudes of data for 
aligned samples are not exact because we cannot exactly esti- 
mate the quantity of the powder contained in the resin. True 
magnitude for H \\ c should be larger than that for H _L c 
because of the c-axis alignment. Inset: A schematic view of 
the crystal structure of (Sr/Pb)Ni2V208- 



The activated temperature dependence of magnetic 
susceptibility is a clear signature of a nonmagnetic 
ground state. The distinct chain-arrangement of mag- 
netic S = 1 Ni 2+ ions in the crystal structure suggests 
that we are in fact dealing with a Haldane antiferro- 
magnet. The behaviour of high-field magnetization is 
also typical of a Haldane-gap system R. The starting 



point for describing the magnetism of PblN^X^Os should 
therefore be a nearest-neighbor Heisenberg Hamiltonian 
with a possible uniaxial single-ion anisotropy (tetragonal 
crystal structure!): H = J £\ S 4 • S m + D ) 2 . For 

not too large D the ground state is a Haldane singlet. 
The lowest energy excitations are a single longitudinal 
mode and a doublet of transverse-polarized Haldane ex- 
citations. Using Eq. (2.14) in Ref. S, relating the criti- 
cal fields i?c and Hc^ to the doublet and singlet gap 
energies Aj_ and A||, for PbNi 2 V20s, we obtain Aj_ = 
2.2 meV and A|| = 1.2 meV. The intrachain exchange 
constant J can be deduced from the high-temperature 
part of the measured x(T) curves §. For PbNi 2 V 2 8 
this analysis yields J = 95 K = 8.2 meV. Already we see 
that the mean gap energy (2Aj_ + Ay)/3 = 1.87 meV is 
substantially smaller that the value 0.41J w 3.36 meV, 
expected for noninteracting chains J|,|| ■ The most likely 
source of this discrepancy is substantial interchain ex- 
change interactions. In their presence the excitation en- 
ergies become dependent on the momentum transfer per- 
pendicular to the chain axis. A_l and A|| in this case 
represent the global minima in the 3D dispersion, and 
cannot be uniquely related to J. 
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FIG. 2. Magnetization curves measured in PbNi2V208 at 
T — 4.2 K for magnetic fields parallel and perpendicular to 
the chain axis. 



Dynamic spin correlations in PbN^X^Os were studied 
in a series of inelastic neutron scattering experiment on 
a 10 g powder sample at the DruchaL 3-axis spectrom- 
eter (SINQ spallation source, PSI, Switzerland). Neu- 
trons of fixed final energy Ef = 8 meV were used with 
a (guide)-80'-80'-(open) collimation setup. An entire se- 
ries of constant--E scans was collected in the energy range 
2-8 meV to map out the low-energy part of the mag- 
netic dynamic structure factor S(Q,u>). Additional data 
were collected in constant-Q scans. A typical data set 
measured at a momentum transfer just exceeding the ID 
antiferromagnetic zone-center = An / c at T = 2 K is 
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shown in the inset in Fig. |3|. This and all other const-Q 
scans clearly show the presence of a spin gap of roughly 
3.5 meV. The main panel in Fig. || shows a constant- .E 
scan measured at hcu = 4 meV, i.e., just above the gap 
energy. In order to analyze these data we have calculated 
the powder average of the 3D dynamic magnetic struc- 
ture factor for an isolated screw-shaped Haldane chain, 
as appropriate for PbNi 2 V 2 0s jnj. The relevant param- 
eters of this cross section are the two gap energies Aj_ 
and A||. Although this model (dashed lines in Fig. |^) 
can reproduce the measured const- Q scans fairly well and 
predicts the asymmetric feature at Q^, it cannot account 
for the sharp peak seen at Q « 37r/c in constant-i? scans. 




Q (units of c*) 

FIG. 3. Typical constant-Q (inset) and constant-_B 
(main panel) inelastic neutron scattering scans measured in 
PbNi2V208 powder at T = 2 K. Lines are model fits, as de- 
scribed in the text. 

A better agreement with the experiment could only 
be obtained by taking into account interchain interac- 
tions. We assumed a weak Heisenberg coupling between 
nearest-neighbor Ni chains and treated it on the RPA 
level. A detailed description of this analysis is beyond 
the scope of this Letter, and will be published elsewhere 
p0| . Compared to the noninteracting-chain model only 
one additional parameter, namely the interchain coupling 
constant J± is introduced. The total number of vari- 
ables is minimized by fixing the lowest excitation ener- 
gies to the gap values deduced from high-field measure- 
ments. An excellent global fit to the data is obtained 
with Ji = 0.096 ± 0.003 meV from least squares refine- 
ment (solid lines in Fig. ||) . From this value and the mea- 
sured critical fields the intrinsic gap energies for nonin- 
teracting chains can be determined: A® =4.1 meV and 
Aii "* = 3.7 meV, respectively. Using the numerical re- 
sults from Ref . , for the intrachain coupling constant J 
and the anisotropy parameter D we obtain J — 9.5 meV 
and D = —0.23 meV (easy-axis), respectively. The for- 
mer value is in reasonably good agreement with the high- 



T susceptibility estimate J = 8.2 meV. 

The results presented above enable us to place the 
new material on the Sakai-Takahashi D — J± phase di- 
agram @. For PbNi 2 V 2 8 we have D/J « -0.024, 
J±/J = 0.01 and z = 4, z being the coordination num- 
ber for interchain interactions. Using Fig. 8 in Ref. [0, 
we find that PbNi 2 V 2 0s is positioned in the Haldane 
phase, just below the borderline with the ordered Ising- 
like phase. In other words, by sheer chance the system is 
almost critical towards the formation of long-range order. 
As will be reported elsewhere unlike PbNi 2 V 2 0s, 
the isostructural SrNi 2 V 2 0s system falls on the oppo- 
site side of the transition line and actually ordered in 3D 
space in a weak-ferromagnetic structure below TV = 7 K. 
A simple substitution of nonmagnetic Sr for nonmagnetic 
Pb slightly alters the coupling geometry, increasing in- 
terchain interactions or single-ion anisotropy beyond the 
critical values. 

Having characterized the pure PbNi 2 V 2 0s compound, 
we now turn to the main subject of this paper, describ- 
ing the emergence of Neel-like magnetic LRO upon in- 
troduction of spin vacancies in this Haldane-gap system. 
Nonmagnetic sites are introduced into the spin chains 
by substituting 5 = Mg 2+ ions for the 5 = 1 Ni 2+ . 
The measured temperature dependences of magnetic sus- 
ceptibility for a series of Pb(Nii_ a; Mg :I ;) 2 V 2 08 powder 
samples with x ranging from to 0.120 are shown in 
Fig. ^. Except for the x = compound, an increase of 
susceptibility is observed at low T in all samples studied. 
The enhancement scales with Mg concentration, and can 
be attributed to the presence of free 5=1/2 spins at 
the ends of the Ni chains severed by the spin vacancies 
p2| . In this scenario the number of liberated spins is 
twice the number of substitutions. To check this scaling 
we compared the temperature dependence of magnetic 
susceptibility in Mg 2+ (5 = 0)- and Cu 2+ (5 = 1/2)- 
substituted PbNi 2 V 2 0s. The replacement of a Ni 2+ ion 
by a Cu 2+ ion also produces two end-chain spins. How- 
ever, these two spins couple antiferromagnetically to the 
Cu-spin, producing a ground state configuration of to- 
tal 5 = 1/2. The total number of added free spins in 
this case is equal to the total number of added 5=1/2 
ions. Experimentally, the measured x(^) curve for 2%- 
Cu-doped PbNi 2 V 2 08 sample was found to be virtually 
indistinguishable from that for a 1%-Mg-doped sample. 

The central result of this work is the observation 
of a sharp decrease in the magnetic susceptibility of 
Pb(Nii_ x Mg x ) 2 V 2 8 below T c » 3.5 K for x > 0.020. 
The slope of x(T) changes abruptly with reversing its 
sign. This behavior is a clear indication of a magnetic 
phase transition. Its nature is revealed in experiments 
on an aligned powder sample with x — 0.030 (Inset 
of Fig. ||). In this case the transition point is about 
T c = 3.0 K. Below this temperature a large anisotropy 
of susceptibility is observed. For H || c, x(T) decreases 
abruptly below T c . In contrast, for H _L c it remains 
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practically constant. This behavior is typical of a Neel 
antiferromagnet with spins ordering along the c direction, 
that, as we have shown previously, is the magnetic easy 
axis of the system. 

The effect of spin-vacancy substitution has been stud- 
ied experimentally (see, for example, Ref. 12 and theo- 
retically 0) for a number of Haldane-gap systems. Apart 
from the appearance of free end-chain spins and the re- 
sulting paramagnetic contribution to the low-T suscep- 
tibility, severing the chains leads to a reduction of the 
dynamic spin correlation length. PbNi2V20s is, to our 
knowledge, the first example of a Haldane AF that orders 
magnetically upon doping. The key of course is in the in- 
terchain interactions: to order magnetically the liberated 
end-chain spins need to be coupled in 3D space. As we 
have demonstrated, interchain coupling is quite substan- 
tial in PbNi2V2C>8. The possibility of vacancy-induced 
LRO in a Haldane-gap system was first qualitatively dis- 
cussed by Shender and Kivelson In particular, the 
authors suggested that there is a critical impurity con- 
centration below which LRO does not occur. As to this 
point more detailed low-T studies are required. 
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FIG. 4. Temperature dependences of magnetic susceptibil- 
ity measured in P^Nii-^Mg^V^Os (0 < x < 0.120) pow- 
der samples. The inset shows the data for a highly aligned 

Pb(Ni .970Mg .030)2V 2 O8. 

The magnetic transition that we observe in 
Pb(Nii_ a Mg a )2V208 should not be confused with that in 
the extensively studied Haldane system CSN1CI3 
Unlike in PbN^X^Os, direct interchain interactions in 
CsNiCl3 are stronger than the critical value, and even 
the undoped material has a magnetically ordered ground 
state. In fact, the introduction of Mg 2+ vacancies leads 
to a decrease of the ordering temperature p5[ , in striking 
contrast with the behavior of PbNi2V20s. Rather, the 
situation in Mg-doped PbN^A^Os is similar to that in 
the spin-Peierls CuGe0 3 , doped with Zn Mg |l7) or 
Ni fill, and the two-leg ladder system SrCu203 doped 



with Zn Although these materials are half-integer- 
spin antiferromagnets, their ground states fall in the 
same universality class as the Haldane state, and they 
both have a gapped magnetic excitation spectrum. Just 
as we propose for PbNi 2 V20g, the onset of long-range 
order is believed to be caused by interchain interactions 
that couple free spins liberated by impurity substitution. 
Doping-induced long-range magnetic order thus appears 
to be a universal feature of quasi-lD singlet quantum 
antiferromagnets with the excitation gap. 

In summary, we have shown that Pbf^X^Os is a 
Haldane-gap system with substantial interchain inter- 
actions and is almost critical towards the formation of 
long-range order. For the first time long-range magnetic 
ordering was induced in a Haldane antiferromagnet by 
spin-vacancy substitution on the 5=1 sites. 

We would like to thank T. Masuda for the coorperation 
in the various aspects of this work. This work was sup- 
ported in part by Grant-in- Aid for COE Research "Phase 
Control of Spin-Chargc-Photon (SCP) Coupled System" 
from the Ministry of Education, Science, Sports, and Cul- 
ture. Work at Brookhaven National Laboratory was car- 
ried out under Contract No. DE-AC02-98CH10886, Di- 
vision of Material Science, U.S. Department of Energy. 



[1] F. D. M. Haldane, Phys. Lett. A93, 464 (1983); Phys. 

Rev. Lett. 50, 1153 (1983). 
[2] T. Sakai and M. Takahashi, Phys. Rev. B 42, 4537 

(1990). 

[3] E. F. Shender and S. A. Kivelson, Phys. Rev. Lett. 66, 
2384 (1991). 

[4] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. 
Rev. Lett. 59, 899 (1987); Commun. Math. Phys. 115, 
477 (1988). 

[5] B. Wanklyn, J. Mater. Sci. 13, 89 (1978). 

[6] R. Wichmann and Hk. Miiller-Buschbaum, Rev. Chim. 

Miner. 23, 1 (1986). 
[7] K. Katsumata et al., Phys. Rev. Lett 63, 86 (1989); Z. 
Honda, H. Asakawa, and K. Katsumata, Phys. Rev. Lett. 
81, 2566 (1998). 
[8] O. Golinelli, Th. Jolicoeur, and R. Lacaze, J. Phys.: 

Cond. Matt. 5, 7847 (1993). 
[9] V. Gadet et at, Phys. Rev. B 44, 705 (1991). 
[10] A. Zheludev et at, to be published. 
[11] Y. Uchiyama et at, to be published. 
[12] M. Hagiwara et al, Phys. Rev. Lett. 65, 3181 (1990). 
[13] W. J. L. Buyers et al, Phys. Rev. Lett. 56, 4 (1986); 

Phys. Rev. B 38, 543 (1988). 
[14] K. Kakurai et al, J. Phys. (Paris) Colloq. 49, C8-1433 
(1988). 

[15] S. S. Sosin et al, JETP 85, 114 (1997). 
[16] M. Hase, I. Terasaki, and K. Uchinokura, Phys. Rev. 
Lett. 70, 3651 (1993); M. Hase et al, Phys. Rev. Lett. 



4 



71, 4059 (1993); Physica B 215, 146 (1995). 
[17] T. Masuda et ai, Phys. Rev. Lett. 80 4566 (1998). 



[18] N. Koide et al, |:ond-mat/9805059| (1998). 

[19] M. Azuma et ai, Phys. Rev. B 55, R8658 (1997). 



5 



